The first jamming crossover: geometric and mechanical features 
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The jamming transition characterizes athermal systems of particles interacting via finite range 
repulsive potentials, and occurs on increasing the density when particles cannot avoid making con- 
tacts with those of their first coordination shell. We have recently shown [M. Pica Ciamarra and 
P. Sollich, arXiv: 1209. 3334 that the same systems are also characterized by a series of jamming 
crossovers. These occur at higher volume fractions as particles are forced to make contact with 
those of subsequent coordination shells. At finite temperature, the crossovers give rise to dynamic 
and thermodynamic density anomalies, including a diffusivity anomaly and a negative thermal ex- 
pansion coefficient. Density anomalies may therefore be related to structural changes occurring at 
the jamming crossovers. Here we elucidate these structural changes, investigating the evolution of 
the structure and of the mechanical properties of a jammed system as its volume fraction varies 
from the jamming transition to and beyond the first jamming crossover. We show that the first 
jamming crossover occurs at a well defined volume fraction, and that it induces a rearrangement of 
the force network causing a softening of the system. It also causes qualitative changes in the normal 
mode density of states and the spatial properties of the normal mode vectors. 

PACS numbers: 61.43.Er; 62.10.+S; 61.20.Ja 



I. INTRODUCTION 



Athermal soft matter systems with particles interact- 
ing via repulsive contact forces, such as granular materi- 
als, foams and emulsions, undergo a percolating transi- 
tion known as the jamming transition when, on increas- 
ing the density, particles can no longer avoid making con- 
tacts. We refer to Refs. [ll, Q for recent reviews. At 
the jamming transition a percolating network of con- 
tact forces appears, and the system acquires mechani- 
cal rigidity. We have recently shown that the jamming 
transition is the first of a series of higher-order jamming 
crossovers 0] . The origin of these crossovers can be un- 
derstood by considering the radial distribution function 
of dense particle systems, which allows one to define a 
sequence of coordination shells. The jamming transition 
marks the point where each particle contacts those of 
its first coordination shell. The jamming crossovers then 
occur at higher density as particles interact with those 
of subsequent coordination shells. Accordingly, at zero 
temperature these crossovers are characterized by an in- 
crease of the mean number of contacts per particle. This 
is accompanied by a change in the mechanical response 
of the system, which becomes more non-affine. At finite 
temperature this leads to density anomalies, such as an 
increasing diffusivity upon isothermal compression, and a 
negative thermal expansion coefficient @ . These anoma- 
lies are analogous to those observed in water and in some 
other network forming liquids 043 ■ Accordingly, the in- 
vestigation of the jamming crossovers may shed light on 
the microscopic origin of these phenomena. In particular. 
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the observed correlation between the jamming crossovers 
and the density anomalies suggest that these anomalies 
might be related to the structural changes occurring at 
the crossovers. 

In this paper we investigate the structural changes oc- 
curring at the first jamming crossover in a bidisperse mix- 
ture of disks interacting via a finite-range, repulsive har- 
monic potential. We show that there is a well defined 
volume fraction (j)'j at which the formation of contacts 
with particles of the first coordination shell ends, and 
a larger well defined volume fraction 0j > 0j at which 
the formation of contacts with particles of the second co- 
ordination shell begins; this is the location of the first 
jamming crossover. We provide evidence that in the vol- 
ume fraction range 0j) two particles are in contact if 
and only if they are also Voronoi neighbors, and the aver- 
age contact number per particle Z is then constant and 
fixed io Z = ^Euicr = 6 by Euler's theorem for planar 
graphs. The formation of contacts between non- Voronoi 
neighbors, which occurs for </) > induces mechanical 
instabilities and avalanches that restructure the whole 
force network, with the net effect of producing a large 
number of weakly compressed bonds. These structural 
changes allow us to rationalize the softening of the sys- 
tem on compression, as well as the volume fraction de- 
pendence of the normal mode density of states and the 
spatial structure of the normal mode vectors. 



II. NUMERICAL MODEL 

We investigated the zero temperature properties of a 
50:50 mixture of soft disks, with diameters D\ and Dg = 
Di/1.4, and equal mass M . Two particles i and ?' interact 
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FIG. 1: (color online) Volume fraction dependence of the 
mean contact number (a), its volume fraction derivative (b), 
and the bulk modulus (c). 

via a harmonic potential when in contact, 

= (1) 

where Sij = Dij — Vij is the overlap between the particles. 
This is expressed in terms of Dij = {Di+Dj)/2, the aver- 
age of the two particle diameters Di and Dj , and the dis- 
tance Tij = |ri — rjl between the particles; 9{x) in Eq. ^ 
is the Heaviside function. We will take M, D\ and k as 
our units of mass, length and energy, respectively. For 
each volume fraction considered, 50 independent jammed 
configurations were generated by minimizing the energy 
via the conjugate-gradient algorithm starting from 
an initial configuration with randomly placed particles. 
We considered values of volume fraction separated by at 
most A(/) = 5 X 10~^, and used measurements taken at 
successive values of the volume fraction to estimate vol- 
ume fraction derivatives @. The results reported below 
concern systems of = 10^ particles unless stated oth- 
erwise. 




FIG. 2: (color online) Dependence oi \Z — ZEuicrl, with 
^Euicr = 6, on the volume fraction (panel a) and on A(j!> (panel 
b), where = — (f> for </) < (j>'j, and = (ji — for 
(j) > Dependence of 9i and 62 (see text) on (j) (panel c) 
and on (panel d). The figure clarifies that Z = ZeuIcf in 
the range <f)'j < (f) < (j>]j, as two particles are in contact if and 
only if they are Voronoi neighbors. 

while Z grows monotonically on compression (Fig. [TJi), 
the rate of formation of new contacts oscillates (Fig. [TId) . 
These oscillations can be interpreted as higher order jam- 
ming crossovers. Indeed, while at the jamming transition 
dZ/d(f) diverges when particles make contacts with other 
particles in their first coordination shell, at higher volume 
fraction dZ/dcf) increases again as particles form contacts 
with subsequent coordination shells. The oscillations of 
dZ/dcf) are correlated with those of the bulk modulus K, 
illustrated in Fig. [TJ;, suggesting a relation between the 
geometric and the mechanical properties of the system. 
In particular, the jamming crossovers appear to induce 
a softening of the system, as K decreases when dZ/d(j) 
increases. The finite temperature counterparts of this 
softening include diffusivity anomalies, and a negative 
thermal expansion coefficient Q. 



III. JAMMING CROSSOVERS 

Fig. [TJ; illustrates the volume fraction dependence of 
the mean contact number per particle Z, its volume frac- 
tion derivative dZ / dcj), and the bulk modulus K. In agree- 
ment with previous results for harmonic disks [H, Q j a-t 
the jamming transition, ~ 0.84, the mean contact 
number jumps to the isostatic value Ziso = 2d = 4, 
and dZ/dcf) diverges. In addition, on approaching the 
jamming transition from above, Z varies as Z — Ziso oc 
{(j) — The bulk modulus, finally, changes discon- 

tinuously at the transition. 

Here we focus on the volume fraction dependence of 
these quantities at higher volume fractions, and note that 



IV. THE FIRST CROSSOVER: GEOMETRIC 
PROPERTIES 

A. The crossover volume fraction, 

In the affinc approximation, the rate of formation of 
new contacts is simply related to the radial distribution 
function. Each jamming crossover would then be iden- 
tified with the volume fraction where the formation of 
contacts with particles of a new shell starts. It is possi- 
ble that the process of formation of such contacts might 
terminate before the next jamming crossover begins. In 
the intervening range of volume fractions, Z is then con- 
stant and dZ/d(j) = 0: all contacts with particles of a 
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given shell are established, and no contacts with parti- 
cles of the subsequent shell form. However, this scenario 
is not very common as successive crossovers usually mix, 
leading to a rate of formation of new contacts that is 
generically nonzero, dZ/dcj) > 0. 

Harmonic disks in two dimensions are peculiar in this 
respect, as the formation of contacts in the first neighbor 
shell, which starts at the jamming transition (jjj, ends 
at a volume fraction (jy'j well below the first jamming 
crossover at (h\ > 6',. We show that this is the case 



by investigating the volume dependence oi \Z — 6\, mo- 
tivated by the fact that Fig. [T^ suggests the existence of 
a volume fraction range where the number of contacts is 
constant at Z = 6. Fig. [5^1 shows that \Z — Q\ vanishes 
in a finite volume fraction range, which extends from (j)'j 



to 



J. We estimate the two volume fractions delineat- 
ing this regime by assuming a power-law dependence of 
\Z — 6\ on A(j). The latter is defined as A0 = j — 
when Z < 6, and as A(j) = (f> — when Z > 6. From 
the numerical fits in Fig. [2t we determine ^ 1.27 and 
<j)j ~ 1.40. We note, however, that accurate estimates 
of these volume fractions are difficult to obtain because 
of the high exponents characterizing the relevant power 
laws. In addition, it is worth remembering that the val- 
ues of 0j and may be slightly protocol dependent, in 
analogy with the location i 
transition [lOll. 




FIG. 3: (color online) (a) Overlap probability distribution 
Pz{S), for increasing volume fractions beyond the first jam- 
ming crossover, as indicated. A negative 5 corresponds to 
the distance between the surfaces of non-interacting parti- 
cles. The distribution is normalized across the positive over- 
laps, so that Jjj Pz{S)dS = Z. (b) To highlight the non-afiine 
rearrangements occurring at the jamming crossovers, we plot 
the actual overlap probability distribution at volume fraction 
(f) — 1.55 and the one resulting from an affine compression of 
the system from ((> = 1.40 to = 1.55 (solid line). 



ij of the conventional jamming 



The value Z ~ Q oi the mean contact number in the 
volume fraction range leading up to the first jamming 
crossover can be rationalized. Indeed, Euler's theorem 
for planar graphs fixes to 6 the average connectivity of 
any tessellation of space in two dimensions, in the limit 
of large system size. Since the simplest tessellation is 
the Voronoi one, we speculate that Z = ^Euier = 6 for 
4>'j < 4> < because in this volume fraction range two 
particles are in contact if and only if they are also Voronoi 
neighbors. To check this hypothesis, we investigate the 
volume fraction dependence of two order parameters re- 
lated to the Voronoi (radical [ll|) tessellation. 



The first one, 9i, is the fraction of Voronoi neighbors 
that are not in contact; the second one, 02, is the fraction 
of contacts that occur between particles which are not 
Voronoi neighbors. Fig. shows that Oi and O2 are 
non-zero for Z < ^Euicr and for Z > Z^uier, and vanish 
at and at respectively. Indeed, their dependence 
on A0 is well described by power laws, as in Fig. [5Jl. 
The results of Fig.[2l3 confirm that, for (f)'j < 4> < 



two 



particles are in contact if and only if they are also Voronoi 
neighbors, and thus rationalize the observed constancy of 
the mean contact number Z in this volume fraction range. 
We note as an aside that this result is not observed for 
all interaction potentials Q, and does not hold in three 
dimensions where Euler's theorem does not constrain the 
mean contact number of the Voronoi tessellation. 



B. Overlap PDF 



We next look at what happens for volume fractions 
above where the mean contact number Z begins 
to rise again with (j). Particles now start to make con- 
tacts with those of their second coordination shell. As 
these new contacts are formed, the entire force network 
undergoes a marked restructuring. To investigate the 
nature of this, we study the volume fraction evolution 
of the overlap probability distribution function, Pz{5). 
We normalize Pz across the interacting particles, so that 
/q^ Pz{5)d5 = Z. K negative overlap corresponds to the 
distance between the surfaces of two particles that are 
not in contact. (This is the reason why we monitor dis- 
tributions of overlaps rather than forces; while for pos- 
itive (5 force and overlap are identical for our harmonic 
potential, forces are zero for all particle pairs with neg- 
ative b and so provide less information on the packing 
geometry.) Figure [3^ shows that the evolution of Pz{S) 
across the first jamming crossover is characterized by the 
growth of a peak at small values of 5 > 0, and by the 
simultaneous decrease of the peak at larger overlaps. 

Such a sharp variation of the overlap probability dis- 
tribution signals the presence of a non-affine evolution of 
the force network on compression. Indeed, in the affine 
approximation the overlap Sij between particles i and j 
evolves as (with the spatial dimension being d = 2 in our 
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case) 



and consequently the overlap distribution around (5 = 
mainly shifts on compression. We illustrate this point 
in Fig. [Sb, where we compare the distributions Pz at 
volume fraction = 1.55 for two different systems. 
One is obtained by energy minimization via our usual 
conjugate-gradient protocol, the other by minimizing 
the energy at = 1.40 and then affinely compressing 
the system to cj) = 1.55. The two distributions differ 
strongly, and the comparison demonstrates that the jam- 
ming crossover induces a non-afRne transformation of the 
force network on compression whose main effect is a large 
increase in the number of small overlaps. 

These small overlaps must mainly result from the re- 
laxation of old compressed bonds. Indeed, since at 
(j) = 1.40 the overlap distribution close to i5 = is small, 
Pz{5 ~ 0) ~ 0, we know that at this volume fraction 
there are very few particles that almost touch and so 
can be brought into contact on compression. The fact 
that the evolution of the overlap distribution on com- 
pression is characterized not only by the growth of a 
peak at small overlaps but also by a decrease of the 
peak at large overlaps, further suggests that some com- 
pressed bonds relax. To see this process at work we in- 
vestigate the evolution of the overlaps in a system pre- 
pared at = 1.40 and then compressed quasistatically, 
with energy minimized at each step in volume fraction, to 
= 1.55. In the scatter plot of Fig. 2] we show the over- 
laps of each pair of particles that is 

in contact before and/or after the compression. The fig- 
ure illustrates that a number of large overlaps at (f) = 1.40 
become small at (j> ~ 1.55. While there also some new 
bonds with small overlaps, the majority of small over- 
laps arises by rearrangement of previously strongly com- 
pressed bonds. This clarifies that the formation of a rel- 
atively small number of contacts with particles in the 
second coordination shell (.Z ~ 8 at </) = 1.55) induces a 
rearrangement of the whole force network, with the latter 
leading to an abundance of small contacts. 

C. Avalanches 

We saw in Fig.|3]that, when a system beyond the first 
jamming crossover is compressed in a quasistastic man- 
ner, particle overlaps change significantly, sometimes by 
amounts of the order of the overlap itself. This suggests 
that these changes do not originate from a continuous 
deformation of the system, but rather from large scale 
catastrophic rearrangements, i.e. avalanches related to 
mechanical instabilities. Indeed, on increasing the vol- 
ume fraction one changes the linear system of equations 
that governs the mechanical equilibrium of the system, 
and may therefore induce a mechanical instability. These 
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FIG. 4: (color online) Scatter plot of the overlaps between 
contacting particles at = 1.40, and at <j} — 1.55. Points 
below the diagonal in the top right quadrant correspond to 
overlaps whose magnitude is decreased on compression. The 
top left quadrant shows newly formed contacts, and the bot- 
tom right quadrant contains contacts that are destroyed. The 
figure clarifies that the peak at small overlaps in Pz{S) at 
(fi = 1.55 is mainly the result of the relaxation of larger over- 
laps that existed at (p = 1.40. 



instabilities are well known to occur in athermal systems 
under quasistatic shear, where they lead to sudden drops 
of the shear stress. To our knowledge avalanches in com- 
pression, which are expected to lead to discontinuities in 
the pressure, have never been reported. 

We have confirmed this avalanche scenario by investi- 
gating the volume fraction dependence of the pressure, 
P, and of the energy, E, of a system of = 10"^ par- 
ticles that is compressed quasistatically from (p ~ 1.40 
to (j> = 1.50 in increments of Scj) = 10"^. Fig. [5] shows 
that the volume fraction dependence of the pressure is 
characterized by a large number of discontinuities; these 
are more clearly visible in the derivative dP/dcj). Dis- 
continuities are also present in the energy as a function 
of (f). They are small and therefore difficult to see by 
eye on the scale of the graph, but can again be detected 
by looking at the volume fraction dependence of dE/dcf). 
Note that all avalanches induce a decrease of the elastic 
energy of the system, as expected. On the other hand, 
even though most avalanches also induce a decrease in 
pressure, there are some that have the opposite effect. 
Since on compressing the system the average value of the 
pressure continues to increase, the compression-induced 
dynamics clearly does not reach a steady state, providing 
an important distinction to the case of quasistatic shear. 
In this respect, avalanches in compression bear more re- 
semblance to those associated with Barkhausen noise in 
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FIG. 5: (color online) Evolution of the pressure (left panels) 
and of the energy (right panels) of a system compressed qua- 
sistatically from cj) = 1.40 to (f) = 1.55. The lower panels show 
the corresponding derivatives, evaluated with volume fraction 
steps of dcf) = 10"^ The fi gure demonstrates the existence of 
avalanches that cause a rearrangement of the force network 
of the system. 



ferromagnetic systems. 



D. Ordering 

Although global ordering is not expected in bidisperse 
systems such as the one studied here, it is still possible 
that compression to high volume fractions may change 
the local order around particles. We next show that 
such a change does indeed occur at the first jamming 
crossover. To do this we measure local ordering via an 
extension of the local bond-orientational order parame- 
ter probing n-fold symmetry in two dimensions |12| . This 
extended parameter is defined as follows: 



1 ^ 1 
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(2) 



Here the first sum runs over all particles z, and the sec- 
ond one over all rii overlaps i5y- of particle i with other 
particles. The angle Qij measures the orientation of each 
bond. When g = 0, all bonds contribute equally to 
which then reduces to the commonly investigated local 
ordering parameter [T^]. When 7^ the weight of a 
bond depends on its magnitude: is mainly affected 
by small overlaps if g < 0, and by large ones if g > 0. 

Considering that the average contact number per par- 
ticle is Z = 4 at the jamming transition, Z = 6 before 
the first jamming crossover, and that it reaches Z = 8 
at (/) ~ 1.55, we have investigated the volume fraction 
dependence of the usual (g = 0) bond order parameter 
for n = 4,6 and 8. Fig. [6] shows that and ^% decrease 
monotonically with the volume fraction, and that this 



FIG. 6: (color online) Volume fraction dependence of the or- 
der parameter for different n and g, as indicated. 



decrease is fastest close to the first jamming crossover, at 
0j ~ 1.40. The order parameter $q, on the other hand, 
is non-monotonic in 0, exhibiting a very small maximum 



close to 



1.40 as well as a more pronounced one at 



larger (\). 

This result is somewhat surprising, as particles with a 
local environment exhibiting an 8-fold symmetry are also 
characterized by a 4-fold symmetry, which one might ex- 
pect to be revealed by Vl/g. To clarify this we investigate 
the influence of small and large overlaps on the order pa- 
rameter, by looking at V^lj ^^^d respectively. The 
results for these parameters in Fig. [HI which show a clear 
peak in ^'l2i demonstrate that at the jamming crossover 
small bonds arrange locally with 4-fold symmetry as in a 
square lattice. There is a much broader peak in VPg, sug- 
gesting that the first jamming crossover also leads to an 
increase of the square symmetry of strongly compressed 
bonds. 

We explicitly illustrate the local square ordering of 
bonds induced by the jamming crossover, focussing on 
a system with volume fraction <^ = 1.55. At this value of 
the volume fraction the overlap distribution has a min- 
imum at (Jinin ~ 0.05, and a maximum at (Smax — 0.24, 
as shown in Fig. [3] This suggests labelling an overlap b 
as small if (S < (5niin, and as large \i b > (5inax- A visu- 
alization of the contact network restricted to such large 
and small overlaps, shown in Fig. [71 clearly reveals the 
presence of local order. One locally identifies two inter- 
penetrating square networks, at a 7r/2 angle, and made 
respectively of bonds with small and of large overlaps. 



E. Softening 

Next we look at the link between the non~affine re- 
arrangement of the force network occurring at the jam- 
ming crossover and the observed softening of the sys- 
tem, particularly the decrease of the bulk modulus. That 
there should be such a link can be seen by considering 
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FIG. 8: (color online) Volume fraction dependence of the av- 
erage value of the overlap between particles of different size. 



FIG. 7; (color online) Partial contact network, for (j) — 1.55. 
Thick black and thin red segments correspond to overlaps 
greater that (5max ~ 0.24 and smaller that Smin ~ 0.05, re- 
spectively. 



a monodisperse system of particles of diameter D in two 
dimensions; the generalization to a discrete mixture or 
a genuinely polydisperse system is straightforward. The 
bulk modulus is given by 



Sc) 



(3) 



where the sum extends over all Nc contacts, and V is the 
volume of the system. Assuming small overlaps, 5c <C D, 
this simplifies to 



K<^^—Z4>{5), 



(4) 



where {5) = X^c '^c is the overlap averaged over 

all contacts. To determine when the system softens, i.e. 
dK/dcf) < 0, we assume that the most important effect of 
the jamming crossover is not the variation of the mean 
number of contacts, but rather the restructuring of the 
force network, and therefore approximate Z as constant. 
With this assumption one finds that 



dK 

^ < 



(5) +3 



d{d) d^{S) 



< 0. 



(5) 



Accordingly, a system softens when the average overlap 
decreases quickly: the negative contribution of the sec- 
ond and possibly the third term in the above expression 
may then overcome the positive contribution of the first 
term. In the bidisperse case, the softening condition is 
analogous to Eq. ([5]), but involves the average overlaps 
between the different particle species, large-large, large- 
small and small-small. Fig. \E\ shows that these average 
overlaps have a non monotonic volume fraction depen- 
dence, and that their maximum values occur roughly at 



the jamming crossover. Accordingly, the rearrangement 
of the force network at the jamming crossover, which 
leads to the formation of a large fraction of small con- 
tacts, explains via Eq. ([5]) the softening of the system. 

We note as an aside that in the above approximation 
the pressure of the system scales as P cx 4>{6). A decrease 
of (5) could therefore lead to an instability of the system if 
d{6)/d(f> < —{5) 14). We do not observe such an instability 
in our system. 



F. Normal modes 

One final perspective on the first jamming crossover is 
via the density of normal modes D{oj) and the spatial 
structure of the normal mode vectors. These contain in- 
formation on the mechanical and the dynamical proper- 
ties of the system. In systems of particles interacting via 
repulsive potentials, D{uj) has been investigated mostly 
close to the jamming transition. At higher densities little 
is known about its variation with </> and about the spatial 
properties of the normal modes. 

Here we describe the evolution of the density of states 
D{uj) in a volume fraction range extending from the jam- 
ming transition to beyond the first jamming crossover. 
This evolution is characterized by three different regimes, 
illustrated in the three panels of Fig.[5J changes occurring 
close to the jamming transition, between the jamming 
transition and the first jamming crossover, and above 
the first jamming crossover. 

The density of states for — 0.85 (Fig. [S^), shows 
features typical of systems close to the jamming tran- 
sition. Debye scaling D{io) cx a; is obeyed up to a 
characteristic frequency w*. At higher frequencies D{uj) 
is first constant, and then grows towards a maximum 
value. If we now increase the volume fraction, w* in- 
creases and the plateau disappears. In addition, a new 
peak emerges, so that D{uj) becomes characterized by 
two peaks. Fig. I^b illustrates that, on further increas- 
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FIG. 9: (color online) Evolution of the density of states with 
volume fraction close to the jamming transition (a), between 
the jamming transition and the first jamming crossover (b), 
and around the jamming crossover (c). 



ing (j), D{uj) mainly evolves via a shift of the first peak 
towards smaller frequencies, while decreasing in height. 
This is in line with the approach to the first jamming 
crossover: the characteristic frequency (of the first peak) 
decreases again, making the system similar to one near 
the conventional jamming transition. Finally, above the 
first jamming crossover, the first peak decreases in height 
as we move away from the first jamming crossover. A new 
peak emerges at a slightly larger frequency, as in Fig. [Ht, 
as a precursor to the approach to the second jamming 
crossover. 

Next we wish to understand the spatial structure of the 
normal mode vectors. Close to the jamming transition, 
conventional parameters such as the participation ratio 
and the bond stretching parameters have not provided 
a clear answer to this question (l3l - [l5l |. Here we show 
that some useful information can be extracted at higher 
volume fractions from an extended bond-stretching pa- 
rameter that is designed to probe the behavior of overlaps 
of a given magnitude. For a generic normal mode a and 
for a typical overlap S, we compute 



Sa{S) = 





(uf - u-)n,, 






uf - uf 





1/2 



(6) 



where A(2;) indicates the Dirac 5-function. Here uf is 
the polarization vector of particle i in mode a, iiy = 
(r; — rj)/rij, and the sums run over all contacts. The 
usual bond stretching parameter is as above, but lacks 
the delta functions and is therefore S independent. The 
contribution to Sa (S) of a bond 5ij of magnitude 6 is close 
to unity when the bond is stretched or compressed, while 
conversely it is close to zero when the bond rotates. Ac- 
cordingly, values of Sa{S) close to one indicate that all 
bonds of typical size 6 stretch or compress, while val- 
ues close to zero tell us than the majority of the bonds 
rotates. In the actual calculation of Sa{S), data are av- 
eraged over a small range of overlaps 6 and over normal 
modes with frequencies in a range around some specified 
w. We re-emphasize that while the conventional bond- 
stretching parameter characterizes the deformation of the 
whole contact network of mode a, Sa{6) characterizes the 
deformation of contacts with overlaps of magnitude S. 

We show in Fig. [TOl graphs of Sa{S) for different val- 
ues of the normal mode frequency, with separate pan- 
els referring to different values of the volume fraction 
as indicated. One notes first that on increasing the fre- 
quency the typical value of Sa{S) increases for all volume 
fractions, indicating that at high frequencies most of the 
bonds stretch. At volume fraction <j> = 0.85, which is 
slightly above the jamming transition, Sa{S) is roughly 6 
independent. This indicates that small and large bonds 
play a similar role in the deformation of the system. This 
observation can be rationalized at least qualitatively from 
the fact that at such relatively moderate volume fractions 
the overlap probability distribution decays exponentially. 
Accordingly, there is no strong separation between large 
and small overlaps, and in addition the number of large 
overlaps is small. 

At (f) = 1.30, in between the jamming transition and 
the first jamming crossover, Sa{S) depends more notice- 
ably on the value of the overlap S. In particular, at low 
frequencies, Sa decreases with increasing S. Accordingly, 
low frequency modes are characterized by a stretching or 
compression of bonds with small overlaps, and by the ro- 
tation of the larger ones. This is in agreement with the 
observation that the energy cost of compressing a bond is 
lower for bonds with small overlaps than with large ones. 
The S dependence of Sa then changes on increasing w, so 
that at high frequencies bonds with large overlaps stretch 
more. 

By the time we reach volume fraction (j) ~ 1.50, above 
the first jamming crossover, the form of Sa has changed 
significantly, particularly at small frequencies. For in- 
stance, at uj = 0.2, Sa{S) is characterized by a high peak, 
and assumes small values both for small and large over- 
laps. The peak occurs at a value of the overlap close to 
the location of the minimum of Pz{S), ^min = 0.05 (see 
Fig. [3]). The fact that Sa assumes small values both at 
small and large values of S may be linked to the fact that 
above the first jamming crossover small and large over- 
laps are spatially correlated, as suggested by Fig. [71 The 
presence of a peak of Sa Sit S = Jmin could therefore be a 
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FIG. 10: (color online) Extended bond stretching parameter 
Sa{5), as a function of the typical overlap, for different fre- 
quencies shown in the legend. Each panel refers to a different 
value of the volume fraction, as indicated. 




FIG. 11: (color online) Typical eigenvector of a A'^ = 10"^ 
particle system at a; = 0.2, and (j> = 1-3 (left) and (j> = 1.5 
(right). 

consequence of the low frequency modes being character- 
ized by relative rotation of clusters of particles. Bonds 
within these clusters would have either quite small or 
rather large overlaps, and bond stretching would occur at 
bonds with intermediate overlaps that connect these clus- 
ters. The possible presence of counter-rotating clusters 
of particles is intriguingly reminiscent of the rigid unit 
modes which have been related to the density anomalies 
of network forming liquids [l6|. 

The changes of the spatial properties of the low fre- 
quency energy modes as we go beyond the first jamming 
crossover in volume fraction can also be appreciated by 
direct visualization. As illustrated in Fig. [Til below the 
crossover low frequency eigenvectors have a large scale 
swirling structure, while above the crossover they look 
more disordered. 



V. CONCLUSIONS 

We have shown that the first jamming crossover of a 
binary mixture of disks interacting via a repulsive har- 
monic potential occurs at a well defined volume frac- 
tion, and that it induces well defined structural changes. 
These changes are related to global rearrangements fol- 
lowing plastic instabilities, and lead to the formation of 
an abundance of small overlaps. In turn, the large num- 
ber of small overlaps is responsible for a decrease in the 
average value of the overlap, and hence for the softening 
of the system. We have characterized the geometric and 
mechanical changes occurring at the crossover using or- 
der parameters designed specifically to probe bonds with 
different overlaps and hence contact forces. In particu- 
lar, the study of the spatial structure of the low frequency 
eigenvectors suggested that they may share some features 
of rigid unit modes [3 ■ Accordingly, these modes should 
play an important role in causing the density anomalies 
of these systems It would therefore be interesting 
to repeat our zero-tempcrature study in simple models 
of water and of some other liquids with density anoma- 
lies [13. 

An intriguing question raised by our results concerns 
the connection between the geometric and mechanical 
properties of particle systems. Indeed, since connectivity 
strongly influences the mechanical properties of a net- 
work [l8|, then close to the jamming transition it has 
been possible to express the mechanical moduli in terms 
of the excess contact number AZ These scaling re- 

lations predict that the bulk and the shear modulus both 
vary monotonically with AZ. Our results clarify that this 
monotonic behavior breaks down at jamming crossovers 
occuring at higher volume fractions, and therefore prove 
that the mechanical properties are not simply related to 
the connectivity. Simple arguments to explain mechani- 
cal properties in terms of geometric ones away from jam- 
ming should therefore be developed; we speculate that 
they would need to take account of the heterogeneity of 
the force network as a key feature of highly compressed 
particle systems. 
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